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Abstract 

Wc investigate the possibility of replacing the cosmological constant with gradual condensation 
of a scalar field produced during the decay of a superheavy dark matter. The advantage of this 
class of models to the ordinary quintessence is that the evolution of the dark energy and the dark 
energy are correlated and cosmological coincidence problem is solved. This model does not need 
a special form for the quitessence potential and even a simple <p theory or an axion like scalar 
is enough to explain the existence of the Dark Energy. We show that the model has an intrinsic 
feedback between energy density of the dark matter and the scalar field such that for a large volume 
of the parameter space the equation of state of the scalar field from very early in the history of the 
Universe is very close to a cosmological constant. Other aspects of this model are consistent with 
recent CMB and LSS observations. 

1 Introduction 

Quintessence models are alternatives to a Cosmological Constant i.e. a non-zero vacuum energy 
density. They are not however flawless. Even in models with tracking solutions the potential of the 
scalar field must somehow be fine-tuned to explain its smallness and its slow variation until today. In 
addition, many of them can not address the coincidence problem i.e. why the density of Dark Matter 
(DM) and Dark Energy (DE) evolve in such a way that they become comparable just after galaxy 
formation. 

Recently a number of authors have proposed interaction between dark matter and quintessence field 
to explain the coincidence. L.P. Chimento et al. based on an earlier work by L.P. Chimento et 
al. j2j and W. Zimdahl et al. j3] suggest an asymptotic scaling law between density of DE and DM. 
In their model due to a dissipative interaction between dark matter and quintessence scalar field <f> q , 
Pdm/Pq — ► cte. where pd m and p q are respectively DM and scalar field density. Assuming this "strong 
coincidence" P, they find the class of potentials V q ((f> q ) such that the equation of state have a solution 
with scaling behavior. Then, using constraints from nucleosynthesis, they find that this category of 
models have w q > —0.7. This value is marginally compatible with WMAP data and far from publicly 
available SN-Ia data which prefers w q ~ — 1. In another version of the same model, W. Zimdahl 
et al. |3j consider a non-static scaling solution, Pdm/Pq °c { a a/o) n - The model with r\ = 1 solves 
the coincidence paradigm but the standard ACDM fits the SN-Ia data better and their best fit has 
Wq ~ —0.7. 

L. Amendola et al. [I] have extensively studied the interaction of quintessence field and dark matter 
in models with tracking solutions and w q > —1. They show that these models are equivalent to 
a Brans-Dicke Lagrangian with power law potential and look like a "Fifth Force". Modification of 
the CMB anisotropy spectra by such interactions is observable and put stringent constraints on their 
parameters. 

D. Comelli et al. [5] study a model in which the effect of interaction between quintessence scalar and 
dark matter appears as time dependence of DM particles mass. This explains the extreme adjustment 



of dark matter and dark energy densities during cosmological evolution. The coupling between two 
fields increases the parameter space for both and reduces by orders of magnitudes the amount of 
fine tuning. In this respect, as we will see below, their model is similar to what we propose in this 
work. However, there are a number of issues that these authors have not addressed. Cosmological 
observations put strict limits on the variation of fundamental parameters including the DM mass. 
In their model the largest amount of variation happens around and after matter domination epoch. 
The mass variation must leave an imprint on the CMB and large structure formation which was not 
observed. 

In addition to the lack of explanation for coincidence in many quintessence models, it is difficult to 
find a scalar field with necessary characteristics in the frame of known particle physics models without 
some fine tuning of the potential [£]. In general, it is assumed that quintessence field is axion with 
high-order, thus non-renormalizable, interactions with the Standard Model particles (or its super- 
symmetric extension) which is highly suppressed at low energies. However, D. Chung et al. [Zj show 
that any supergravity induced interaction between <p q and other scalars with VEV of the order of 
Plank mass can increase the very tiny mass of the (j) q (m q ~ Hq ~ 10~ 33 eU) expected in many models, 
unless a discrete global symmetry prevents their contribution to the mass. 

In a very recent work, G.R. Farrar & P.J. Peebles jHj study models with a Yukawa interaction between 
DM and quintessence scalar field. Like D. Comelli et aZ.model, this interaction affects the mass of the 
dark matter particles. The general behavior of these models is close to ACDM with some differences 
which can distinguish them. One of the special cases with a 2-component CDM imitates the ACDM 
very closely. Many aspects of this model is similar to the model studied in the present work but 
without considering the source of the intimate relation between DM and DE in contrast (we believe) 
to the present work. Moreover, the necessity of having a very special self-interaction potential for the 
quintessence field is not removed. 

What we propose here is a model for dark energy somehow different from previous quintessence models 
(A preliminary investigation of this model has been presented in [5] ) . We assume that DE is the result 
of the condensation of a scalar field produced during very slow decay of a massive particle. In most of 
quintessence models the scalar field is produced during inflation or reheating period in large amount 
such that to control its contribution to the total energy of the Universe, its potential must be a 
negative exponential (in most cases sum of two exponentials) or a negative power function jHj. We 
show that in the present model very small production rate of the scalar field replaces the fine tuning 
of the potential and practically any scalar field even without a self-interaction has a tracking solution 
for a large part of its parameter space. 

The main motivation for this class of models is the possibility of a top-down solution ^0] HD H2| 
for the mystery of Ultra High Energy Cosmic Rays (UHECRs) HH [15] • If a very small part 

of the decay remnants which make the primaries of UHECRs is composed of a scalar field <f) q , its 
condensation can have all the characteristics of a quintessence field. We show that in this model the 
most natural equation of state for the quintessence scalar is very close to a cosmological constant, at 
least until the age of the Universe is much smaller than the life-time of the Superheavy Dark Matter 
(SDM, WIMPZILLA) which is the origin of the quintessence field. 

Another motivation is the fact that a dark energy with w q < — 1 fits the SN-Ia data better than a 
cosmological constant JE] Q2| ■ Although the sensitivity of CMB data to the equation of state of 
the dark energy is much less than SNs, with 95% confidence WMAP data gives the range — 1±0.22 for 
the w q |19j |2flj . Estimation from galaxy clusters evolution is also in agreement with this range |21j . 
On the other hand, it has been demonstrated that the cosmological equation of state for a decaying 
dark matter in presence of a cosmological constant is similar to a quintessence with w q < — 1 |18j . 
Both observations therefore seem to encourage a top-down solution which explains simultaneously the 
dark energy and the UHECRs. 

Like other models with interaction between DM and DE, the coincidence in this model is solved 



without fine-tuning. Parameters can be changed by many orders of magnitude without destroying the 
general behavior of the equation of state or the extreme relation between the energy density of dark 
energy and the total energy density in the early Universe. 

In Sec[2]we solve the evolution equations for dark matter and dark energy. For two asymptotic regime 
we find analytical solutions for the evolution of <j) q . In Seed] we present the results of numerical 
solution of the evolution equations including the baryonic matter and we show that both approaches 
lead essentially to the same conclusion. We study also the extent of the parameter space. The effect 
of DM anisotropy on the energy density of the dark energy is studied in Sec|lJ We show that the 
perturbation of dark energy in this model is very small and very far from the resolution of present or 
near future observations. The late time decoherence of the scalar field is discussed in Sec l2.ll We give 
a qualitative estimation of the necessary conditions and leave a proper investigation of this issue as 
well as the possible candidates for cj) q to future works. 



2 Cosmological Evolution of a Decaying Dark Matter and a Quintessence 
Field 

Consider that at very early epoch in the history of the Universe, just after inflation, the cosmological 
"soup" consists of 2 species: a superheavy dark matter (SDM) - X particles - decoupled from the rest 
of the "soup" since very early time and a second component which we don't consider in detail. The 
only constraint we need is that it must consist of light species including baryons, neutrinos, photons, 
and light dark matter (by light we mean with respect to X). For simplicity we assume that X is a 
scalar field <f> x . Considering <f> x to be a spinor or vector does not change the general conclusions of 
this work. We also assume that 4> x is quasi-stable i.e. its lifetime is much longer than the present age 
of the Universe. A very small part of its decay remnants is considered to be a scalar field (j) q with 
negligibly weak interaction with other fields. 



The effective Lagrangian can be written as: 
C = I d A x\/—g 



\g^d^ x d v 4> x + ^a^d^g - v{<t> x , </> q , J) 



The field J presents collectively other fields. The term V(4> x , 4> q , J) includes all interactions including 
self-interaction potential for <f) x and (fi q : 

V{<j> x , 4> q , J) = v q {<t> q ) + V x {<j> x ) + g4> x m 4> q n + W{<j> x , cj> q , J) (2) 

The term g4> x m (j)q n is important because it is responsible for annihilation of X and back reaction of 
quintessence field by reproducing them. W(<j) x , 4> q , J) presents interactions which contribute to the 
decay of X to light fields and to cp q (in addition to what is shown explicitly in ©). The very long 
lifetime of X constrains this term and g. They must be strongly suppressed. For n = 2 and m = 2 the 
g term contributes to the mass of 4> x and 4> q . Because of the huge mass of 4> x (which must come from 
another coupling) and its very small occupation number < (f> x 2 >~ 2p x /m x 2 , for sufficiently small g 
the effect of this term on the mass of the SDM is very small. We discuss the role of this term in detail 
later. If the interaction of other fields with <p q is only through the exchange of X (for instance due to 
a conserved symmetry shared by both X and (f) q ), the huge mass of X suppresses the interaction and 
therefore the modification of their mass. This solves the problem of "Fifth Force" in the dark j3] and 
the SM sectors. 

In a homogeneous universe the evolution equations for <f> q and <f> x are: 

dV 

<j) q + 3H0 q + — = (3) 
dV 

(j)x + 3H( p x + — = (4) 



where dot means the comoving time derivative. In the rest of this work we treat (j) x and J as classical 
particles and deal only with their density and equation of state. We assume that X particles are 
non-relativistic (i.e. part of the CDM) with negligible self-interaction i.e. 

Vx{<t>x) = \rn x 2 (p x 2 (5) 



Under this assumptions 4> x can be replaced by: 
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(6) 



If X is a spinor, the lowest order (Yukawa) interaction term in <^Q) is g^gipip. In the classical treatment 
of X: 

tfi/j ~ —2- (7) 
m x 

The same argument about the negligible effect of the interaction on the mass of DM and SM particles 
is applied. For simplicity we consider only the scalar case. 

For potential V q (<p q ) we consider a simple (f) A model: 

VM = lm q % 2 + jcj> q 4 (8) 

Conservation of energy-momentum, Einstein and dynamic equations, give following system of equa- 
tions for the fields: 

4> q [<P q + ZH<j> q + m 2 <j) q + X(/) q 3 } = -2g0 q( j) q ( ) + T qPx (9) 

p x + 3Hp x = -{T q + Tj)p x -^g 2 (^-^\ (10) 
pj + 3H(pj + Pj) = Tjp x (11) 

H 2 = {-) =—(p x + PJ + P q ) (12) 
P q = ^rn q 2 <j) q + \n 2 4> 2 + j(f) g 4 (13) 

where (|10j) is the Boltzmann equation for X particles. We calculate its right hand side in the appendix. 
p q ' is the density of quintessence particles (not the classical field <j) q ) with an average energy larger than 
m x in the local inertial frame. Only interaction between these particles contribute to the reproduction 
of SDM. T q and Tj are respectively the decay width of X to (j) q and to other species. In @ we have 
replaces <j> x with its classical approximation from (jSJ). The effect of decay Lagrangian W(4> x , 4> q , J) 
appears as (r 9 + Vj)p x which is the decay rate of X particles (see equation (jl3*|) in the Appendix). 

At very high temperatures when p x 3> Tr' i g 2 m x 3 T, the annihilation and reproduction terms in ()10j) 
are dominant. X particles however are non-relativistic up to temperatures close to their rest mass. 
Quintessence scalar particles at this time are relativistic and therefore their density falls faster than 
SDM density by a factor of a(t). The probability of annihilation also decreases very rapidly. Con- 
sequently, from very early time only the decay term in (|10|) is important. The dominance of an- 
nihilation/reproduction can happen only if the production temperature of X particles i.e. preheat- 
ing/reheating temperature is very high. Such scenarios however can make dangerous amount of grav- 
itinos For this reason, presumably the reheating temperature must be much smaller than m x and 
annihilation dominance never happens. This can not put the production of SDM in danger because it 
has been shown [23] that even with a very low reheating temperature they can be produced. It seems 
therefore reasonable to study the evolution of the fields only when the annihilation/reproduction is 
negligible. Another reason for this simplification is that we are interested in the decohered modes of 



<ft q . When the self-annihilation of X particles is the dominant source of (ft q most of particles are highly 
relativistic and their self interaction doesn't have time to make long-wavelength modes. This claim 
needs however a detail investigation of the process of decoherence which we leave for another work. 

The system of equations (|9|)- (|13j) is highly non- linear and an analytical solution can not be found 
easily. There are however two asymptotic regimes which permit an approximate analytical treatment. 
The first one happens at very early time just after the production of X (presumably after preheating 
jllj |12j ) and the decoherence of (ft q s long wavelength modes. In this epoch <ft q ~ and can be 
neglected. The other regime is when comoving time variation of (ft q is very slow and one can neglect 
4> q . We show that the first regime leads to a saturation (tracking) solution where 4> q — > cte. It then 
can be treated as the initial condition for the second regime when <ft q changes slowly. 

The effect of last term in right hand side of ()10|) as we argued is negligible. The solution of ()10|) is 
then straightforward: 

-r(t-io) ( 



P*(t) = Px (to)t f 1 (t - t0, [^f) (14) 

where r = T q + Tj is the total decay width of X. We consider to to be the time after production and 
decoupling of X. These two times can be very different, but with an extremely long lifetime for X 
and its weak interaction with other species, it is not important which one of them is selected as to- 

After inserting the solution (|14j) and neglecting all the terms proportional to (ft q , equation @ simplifies 
to: 

4> q j t (a 3 4> q ) = T q a 3 (to)pMe- r{t - to) (15) 

and can be solved: 



-cft q (t) = K q (t) - 



2^"'-"^-' \ a {t) 



K q {to) + T qPx (to) f dt^-e-^ 1 -^ 
Jt a{t Q ) 



(16) 



For a oc t k the integral term in (|16|) decreases with time (i.e. (ft q < 0). This means that after a 
relatively short time (ft q is saturated and its density does not change, in other words it behaves like a 
cosmological constant. The numerical simulation in the next section confirms this result. If <j) q was 
a classical field the natural choice for the initial value of the kinetic energy K q (to) was K q (to) = 
assuming a very rapid production of X. However, in reality 4> q is a quantum field and it gets time to 
decohere and to settle as a classical field. The initial value of K q {to) can therefore be non-zero. Its 
exact value can only be determined by investigating the process of decoherence. In any case with the 
expansion of the Universe, its effect on </> g decreases very rapidly because of a~®{i) factor in (|16|) . 

Next we consider the regime where (f) q changes very slowly and we can neglect (ft q and higher orders 
of 4> q . Equation Q gets the following simplified form: 

<ft q (m q 2 cft q + X(ft q 3 ) = -2gcft\(ft q (^A + T qPx (17) 

We expect that self-interaction of (ft q be much stronger than its coupling to X. Neglecting the first 
term in the right hand side of (|17|) . its ^-dependent part can be integrated: 

K^v+i*» 4 )=fw= r ^ (i8) 

which then is easily solved: 

V q (4> q ) = V q (4M) +r ^( t o)/, < ^(^) 3 e- r( ^ ) (19) 

Here V q is the potential energy of cft q . From (|18|) and (|19|) it is clear that the final value of the potential 
and therefore (ft q energy density is driven by the decay term and not the self-interaction. Therefore 
the only vital condition for this model is the existence of a long life SDM and not the potential of (ft q . 



In (|19l) the initial values t' and 4> q (t' Q ) are different from equation (|16|) . They correspond to the time 
and to the value of <j) q in the first regime when it approaches to saturation. Similar to H16|) . the time 
dependence of 4> q in (|19|) vanishes exponentially and the behavior of 4> q approaches to a cosmological 
constant. 

To estimate the asymptotic value of 4> q we assume that a(t) cx t k . Using (|19|) with the additional 
assumption that t s — t' <C 1/T, (t s is the saturation time), we find: 

/ ,/ \ (3fc-l)\ 

-(t) )■ < 2 °> 

V((f>q(t' )) has 90% of its final value, for 
1/2 and: 

(21) 
(22) 

The interesting conclusion one can make from (|2()|1 is that the initial density of SDM, its production 
time, and its decay rate to 4> q which are apparently independent quantities determine together the final 
value of the dark energy density. The long lifetime of SDM is expected to be due to a symmetry which 
is broken only by non-renormalizable high order weak coupling operators. They become important 
only at very large energy scales. These conditions are exactly what is needed to have a small dark 
energy density according to (|20[) . In Sec 01 we see that numerical calculation confirm these results. 

We can also observe here the main difference between this model and other quintessence models. If 
4> q is produced during e.g. the decay of inflaton or from the decay of a short live particle in the early 
Universe, its final density should be much larger than observed dark energy unless either its production 
width was fine-tuned to unnaturally small values or its self-interaction was exponentially suppressed 
with some fine tuning of its rate. 



v(cf> q )-v(<f>M)-^^(i 
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If we define the saturation time as the time when V(<j) q ) 
t s *C t eq with t eq the matter-radiation equilibrium time, k 

t 8 ~ 100*0 

For t s > t e? , fc = 2/3 and: 
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o 



2.1 Decoherence 

Decoherence of scalar fields has been mainly studied in the context of phase transition [23] in a 
thermal system. Examples are phase transition in condense matter [23] j2£], and before, during and 
after inflation in the early Universe [2B] |27| . In the latter case the aim is studying the inflation 
itself, production of defects and the reheating. Decoherence is the result of self-interaction as well 
as interaction between a field (regarded as order parameter after decoherence) and other fields in 
the environment. Long wavelength modes behave like a classical field i.e. don't show "particle-like" 
behavior if quantum correlation between modes are negligible. More technically this happens when 
the density matrix for these modes is approximately diagonal. It has been shown that interaction 
with higher modes is enough to decohere long wavelength modes (see Calzetta et al. |27j for a review). 
The classical order parameter corresponds to these modes after their decoherence. One can consider a 
cut-off in the mode space which separate the system (i.e. long wavelength modes) from environment 
(short wavelengths). The cutoff can be considered as an evolving scale which determines at each 
cosmological epoch the decoherent/coherent modes |27j . 

It has been shown that the decoherence time in a thermal phase transition is shorter than 

the spinodal time i.e. the time after beginning of the phase transition when the scalar field or more 
precisely < cp 2 > settles at the minimum of the potential. The decoherence time in presence of external 
fields (with couplings of the same order as self-interaction) is 



1 

td ~ — 
m 



(23) 



By replacing Minkovski time with conformal time and considering a time dependent cut-off |2f)] E3 
|25|one can show that modes with: 

- 2 + m 2 < H 2 (24) 
gt 

decohere and behave like a classical scalar field. The effect of coupling constant is logarithmic and 
less important. 

If the SDM exists, it is produced during preheating just after the end of the inflation presumably 
at T ~ 10 14 eU — 10 16 eU which correspond to: 

H ~ KT 6 eU - l(T 4 ey (25) 

From (|21|) this time range permits scalars with mass m < 10 -6 eV to decohere. When the size of the 
Universe get larger, 4> q stops decohering. This also helps having a very small dark energy density. If 
the preheating/reheating had happened when the Hubble constant was smaller, them m q also must 
be smaller to have long wavelength modes which can decohere. We will see in the next section that 
in this case the main term in V((f) q ) potential is the self-interaction. Moreover, A can be larger which 
helps a faster decoherence of long wavelength modes. 

The argument given here is evidently very qualitative and needs much deeper investigation. In the 
present work we take the possibility of decoherence as granted and study the evolution of 4> q as a 
classical scalar field. 



3 Numerical Solution 

To have a better understanding of the behavior and the parameter space of this model, we have solved 
equations © to (fT3|) numerically. We have also added the interaction between various species of the 
Standard Model particles to the simulation to be closer to real cosmological evolution and to obtain 
the equation of state of the remnants. This is specially important for constraining the lifetime of 
SDM |28j . Without considering the interaction between high energy remnants and the rest of the SM 
particles specially the CMB, the lifetime of SDM must be orders of magnitude larger than present age 
of the Universe. 



Details of interaction simulation are discussed in (22] and we don't repeat them here. The Boltzmann 
equation for SM species (equation (1) in |28j ^ replaces Because of numerical limitations we switch 
on interactions only from z = 10 9 downward. For the same reason we had to begin the simulation 
of X decay from z ~ 10 14 which is equivalent to a temperature of T = W n eV. The expected 
reheating temperature is model dependent and varies from ~ W 22 eV to ~ W 7 eV. For the time being 
no observational constraint on this large range is available. The change in the initial temperature 
however does not modify the results of the simulation significantly if f q = T q /T is rescaled inversely 
proportional to redshift and to the total decay width T, and proportional to m x . In other words two 
models lead to very similar results for the quintessence field if: 



f q z'T'wii 



f q zTm' x 



(26) 



For the lifetime of X we use the results from |28j and [IB] which show that a lifetime r = 5tq — 50tq 
(to the present age of the Universe) can explain the observed flux of UHECRs as well as cosmic 
equation of state with w q < — 1. In the following we consider r = 5tq. Our test shows that increasing 
r to 50to does not significantly modifies the extent of the admissible parameter space or other main 
characteristics of the dark energy model. We consider only the models with m = 2, n = 2 and 
g = 10~ 15 in ©. The results for 10~ 20 ^ g ^ 10~ 5 are roughly the same as what we present in 
this section and therefore they are not shown. The discussion in SecEJas well as Fig[2]show that the 
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Figure 1: Evolution of quintessence field (left), its derivative (center) and its total energy density 
(right) for Tq = T q /T = 10 -16 (magenta) (see text for details), 5Fo (cyan), lOFo (blue), 50Fo (green), 
lOOro (red). Dash line is the observed value of the dark energy. m q = 10~ 6 eV, A = 10~ 20 . 



contribution of the interaction with the SDM in the total energy density of 4> q is much smaller than 
other terms. 

FigQ shows the evolution of <f) q , its time derivative and its total energy density from the end of X 
production to saturation redshift z s . Here we have used as z s the redshift after which up to simulation 
precision the total energy density of <p q does not change anymore. The result of the simulation is 
quite consistent with the approximate solutions discussed in SecJ^J The final density energy of 4> q 
is practically proportional to T q /T. The latter quantity encompasses 3 important parameters of the 
model: The fraction of energy of the remnants which changes to <j) q , the fraction of energy in the long 
wavelength modes which can decohere and the coupling of these modes to the environment which 
contributes to (fr q yield and to the effective formation redshift of the classical quintessence field <j) q . 
Therefore the effective volume of the parameter space presented by this simulation is much larger and 
the fine-tuning of parameters are much less than what is expected from just one parameter. 
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Figure 2: Evolution of the contribution to the total energy density of cp q for Tq = T q /T = 10~ 16 and 
: Left, m q = 10~ 8 eV and A = 10~ 20 ; Center, m q = 10~ 6 eV and A = 1(T 20 ; Right, m q = W~ 6 eV and 
A = 10~ 10 . Curves are: mass (red), self-interaction (green), kinetic energy (cyan) and interaction with 
SDM (blue). 



FigED shows the evolution in the contribution of different terms of the Lagrangian to the total 
energy of (ft q . Very soon after beginning of production of quintessence field the potential takes over the 
kinetic energy and the latter begins to decrease. The relative contribution of each term and their time 
of dominance as this figure demonstrates, depends on the parameters specially m q and A. Another 
conclusion from this plot is that changing these parameters by orders of magnitude does not change 
the general behavior of the model significantly and for a large part of the parameter space the final 
density of quintessence energy is close to the observed value. This can also be seen in FigJSland Fig0] 
where the evolution of quintessence energy is shown for various combination of parameters. 
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Figure 3: Left: Evolution of total density with redshift for Tq = F q /T = 10 16 (magenta) (see text 
for details), 5Tq (cyan), IOFq (blue), 50Fo (green), IOOFq (red). Dash line is the observed value of 



the dark energy. 
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Right: Relative density of dark energy and CDM as a 



function of T q /T. The x-axis is normalized to Tq = T q /T = 10 



-16 
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Figure 4: Quintessence energy density for: Left, m q = 10 3 eV (cyan), m q = 10 5 eV (magenta), 
m q = 10~ 6 eV (red) and m q = W~ 8 eV (green), A = 10~ 20 ; Right, A = 10~ 10 (cyan), A = 10" 15 , 
A = 10 -20 and A = 10 -25 (green), m q = W~ e eV. The difference between quintessence density for the 
last 3 values of A is smaller than the resolution of the plot. Dash line is the observed energy density 
of the dark energy. 



4 Perturbations 

Large and medium scale observations show that the dark energy is quite smooth and uncorrelated 
from the clumpy dark matter [22]. If DE origin is the decay of the dark matter, the question arises 
whether it clumps around dark matter halos or has a large scale perturbation which is not observed in 
the present data. In this section we investigate the evolution of spatial perturbations in (J) q and show 
that they decrease with time. Another interest in doing such exercise is to investigate any imprint of 
the model on the power spectrum of matter and the CMB anisotropy. 

We use the synchronous gauge metric: 

ds 2 = dt 2 - a 2 (t)(^- - h lj )dx i dx j (27) 

For small spatial fluctuations <j) q (x, t) = 4> q (t) + 5(j) q (x, t) where from now on barred quantities are the 
homogeneous component of the field depending only on t. We define the same decomposition for other 
fields. 

We consider only scalar metric fluctuations h = 5 tJ hij and neglect vector and tensor components. The 



Einstein equation gives following equation for the evolution of h: 

4nG{4$ q 6'<t> q - 25V(<j> q , Px ) + Sp x + Spj + 3SPj) 



1 ■■ a ■ 

-h+ -h 

2 a 



(28) 



where Sp x is the fluctuation of X particles density, 5pj and 5Pj are respectively the collective density 
and pressure fluctuation of other fields. From the Lagrangian the dynamic equation of (p q is: 



hW^sTdutq) + V=gv'(<p q , <p x , J) = o 

This equation and the energy momentum conservation determine the evolution of 5<f> q (x,t): 



(29) 
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Like homogeneous case, we assume that SDM behaves like a pressure-less fluid: 



oo 



p x + 8p x T x 0z = p x 5u 



T x ij = O(5 2 )^0 



r q {8 Px - -^m®) 

(31) 



(32) 



where bu x is the velocity of SDM fluctuations with respect to homogeneous Hubble flow. Interaction 
terms are explicitly included in the energy-momentum conservation equation: 



Px J 2 



4 2 



35p x 2p q '5p q p q 8p 3 



m. 



m q 3 p x 



m q 3 p x 2 



(33) 



The effect of interactions in the right hand side of ()33|) is however very small, first because X particles 
mass is very large and then because only high energy cp q particles contribute to this term and their 
energy decreases with expansion of the Universe much faster than the SDM. The evolution of matter 
fluctuations is then practically the same as the standard ACDM case. 



Using the conservation relation for other components of the energy-momentum in the limit when 
4> q — ► 0, we find the following relation between spatial fluctuation of 54> q and 5u x : 

-V\f q , p x W(8<t> q ) = T q p x 5u x l (34) 

Equation (|3U|) has a meaningful limit when (f> q — > only if 5(p q — > 0. On the other hand, (|34|) shows 
that the divergence of quintessence field fluctuations d l 5(f) q follows the velocity dispersion of the dark 
matter with opposite direction. Their amplitude however is largely reduced due to the very small 
decay width T q . In addition, with the expansion of the Universe, V'((fr q ,p x ) varies only very slightly - 
just the interaction between SDM and <p q will change. In contrast, p x decreases by a factor of a~ 3 (t) 
and even gradual increase of the dark matter clumping and therefore the velocity dispersion bu x |29| 
can not eliminate the effect of decreasing density. We conclude that the spatial variation of <j> q is very 
small from the beginning and is practically unobservable. 



5 Closing Remarks 



Since the original works on the production of superheavy particles after inflation |12j . a number of 
investigations |23j have demonstrated that even with a reheating temperature as low as few MeV the 
production of superheavy particles is possible. We don't discuss here the particle physics candidates 
for <j) q , but for the sake of completeness we just mention that axion like particles are needed or at 
least can exist in large number of particle physics models (see [SJ f° r some examples). The fact 



that <f)q does not need to have very special potential is one of the advantages of this model with 
respect to others and opens the way to a larger number of particle physics models as candidate for 
the quintessence field. 

One of the arguments which is usually raised in the literature against a decaying dark matter is the 
observational constraints on the high energy gamma-ray and neutrino background. In |28j it has been 
shown that if m x > 10 22 eV and its lifetime r > 5to, and if simulations correctly take into account the 
energy dissipation of the high energy remnants, present observational limits are larger than expected 
flux from a decaying UHDM. Consequently the model is consistent with the available data. 

The same fact is applied to the CMB and its anisotropy. The expected CMB distortion is of order 
10 -8 , much smaller than sensitivity of present and near future measurements. As for the expected 
anisotropy in the arrival direction of UHECRs, the data is yet too scarce to give any conclusive answer. 
In the next few years the Auger Observatory [HOI is able to test top-down models for UHECRs which 
is one of the principal motivation for the quintessence model proposed here. 

Although the limit on the amount of the hot DM can not constrain this model, a better understanding 
of its contribution to the total density and its content can help to understand the physics and the 
nature of SDM if it exists. 

Evidently the observation of w q and its cosmological evolution is crucial for any model of dark energy. 
Observation of small anisotropy in the DE density and its correlation with matter anisotropy also can 
be used as signature of relation/inter action between DM and DE. 

For the range of expected masses for (fi q in this model, the high energy component of quintessence field 
is yet relativistic. As we have discussed in Sec[2J the production of this component from annihilation 
has been stopped very early in the history of the Universe and the contribution from decay of X is 
much smaller than the limits on the amount of Hot Dark Matter (as it has been shown in [2B] for hot 
SM remnants). The small coupling of <j) q with SM particles also suppresses the probability of its direct 
detection. However, the detection of an axion-like particle e.g. the QCD axion can be a positive sign 
for the possibility of existence of c^-like particles in the Nature. 

Appendix 

Here we calculate the right hand side of the Boltzmann equation at lowest order of g coupling constant 
for annihilation and reproduction of X particles. 

The Boltzmann equation for X particles is the following: 



P' 



»d^ x \x,p) - r^V^ = L[f] = -(A(x,p) + B(x,p))f( x \x,p) +C(x,p) (35) 

A(x, P ) = ryv (36) 

B(x,p) = -r-4— Y\ f dp x f^(x,p x )Aa xi i = X,<f> q (37) 

C(x,p) = J d Pq dp q f^\x,p q )f^Hx,p' q )A da ^ + ^ X+X (38) 

A = \/ (pi-P2) 2 - m 1 2 m 2 2 (39) 



The function f( x \x,p) is the distribution of X particles. The terms A, B and C are respectively the 
decay, the annihilation (self or in interaction with other species) and production rates. We assume 
that the interaction of X with other fields except (p q is negligible. According to Lagrangian Q with 
n = 2 and m = 2 the lowest Feynman diagrams contributing to annihilation and production are: 




The S matrix for these diagrams is very simple: 

-z(2tt)V (4) (£^) 



and the differential cross-section: 

da _ (2vr) 1 V<5 (4) (E i P, 



=dp 3 dp4, dpi = 



d 3 pi 



(2vr)3^o 



(40) 



(41) 



16\/ (PI-P2) 2 - mi 2 m 2 2 

where gi is the number of internal degrees of freedom. Here we assume that g x = g q = 1. Using the 
relation: 

J p^ 1 . ..p^f(x,p)dp = JpP 1 . ..p^L[f](x,p)dp (42) 
and the definition of energy- momentum tensor T^ v and number density of particles n^, one obtains: 
T ^u v = _ TT ^ Uv _^ g 2f^ j d p 2 /»( x , p2 )_y dp 1 dp 2 pSf^(x,p 1 )f^(x lP2 )8(j )1 +p 2 -2m x ) ) j 

(43) 

Both and have a large peak around the energies close to the mass of X. Therefore: 

'dpf®(*,P)*— < = ^,9 (44) 



■/??.; 



In the case of 4> q the density n q is only the density of particles with an average energy larger than 
m x . Finally from (|43|) one can obtain the evolution equation of p x in a homogeneous cosmology i.e. 
equation (fl7)|) in SecEJ 
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